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Abstract. We study the effect of non-radial motions on 
the mass function, the velocity dispersion function (here- 
after VDF) and on the shape of clusters of galaxies using 
the model introduced in Del Popolo & Gambera (1998a, b; 
1999). The mass function of clusters, obtained using the 
quoted model, is compared with the statistical data by 
Bahcall & Cen (1992a,b) and Girardi et al. (1998), while 
the VDF is compared with the Center for Astrophysics 
(hereafter CfA) data by Zabludoff et al. (1993) for local 
clusters and those of Mazure et al. (1996) and Fadda et 
al. (1996). In both cases the model predictions are in good 
agreement with the observational data showing once more 
how non-radial motions can reduce many of the discrep- 
ancies between Cold Dark Matter (hereafter CDM) model 
predictions and observational data. Finally we study the 
effect of non-radial motions on the intrinsic shape of clus- 
ters of galaxies showing that non-radial motions produce 
clusters less elongated with respect to CDM model in 
agreement with de Theije et al. (1995, 1997) results. 

Key words: cosmology: theory-large scale structure of 
Universe - galaxies: formation 



1. Introduction 

At its appearance the CDM model contributed to ob- 
tain a better understanding of the origin and evolution 
of the large scale structure in the Universe (White et al. 
1987; Frenk et al. 1988; Efstathiou 1990; Ostriker 1993). 
The principal assumptions of the standard CDM (SCDM) 
model (see also Liddle & Lyth 1993) are: 

— a flat Universe dominated by weakly interacting el- 
ementary particles having low velocity dispersion at 
early times. The barionic content is determined by the 
standard big bang nucleosynthesis model (Kernan & 
Sarkar 1996; Steigman 1996; Olive 1997; Dolgov 1997); 

— critical matter density; 

— expansion rate given hy h = 0.5; 
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— a scale invariant and adiabatic spectrum with a spec- 
tral index, n = 1; 

— the condition required by observations, that the fluctu- 
ations in galaxy distribution, {Sp/p)g , are larger than 
the fluctuations in the mass distribution, (Sp/p)p by a 
factor b > 1. 

If this last assumption is not introduced, the pairwise ve- 
locity dispersion is larger then that deduced from observa- 
tions and the galaxy correlation function is steeper than 
that observed (Davis et al. 1985). After the great success 
of the model in the 80's, a closer inspection of the model 
has shown a series of deficiencies, namely: 

— the strong clustering of rich clusters of galaxies, 
^cc('') — (r/25/i~^Mpc)~^, far in excess of CDM pre- 
dictions (Bahcall & Soneira 1983); 

— the overproduction of clusters abundance. Clusters 
abundance is a useful test for models of galaxy forma- 
tion. This is connected to three relevant parameters: 
the mass function, the VDF and the temperature func- 
tion. Using N-body simulations, Jing et al. (1994) stud- 
ied the mass function of rich clusters at 2 = for the 
CDM model concluding, if the density spectrum is nor- 
malized to the Cosmic Background Explorer (hereafter 
COBE) (Smoot et al. 1992) quadrupole Qcobe = 
6.0x 10~®, that the mass function is higher than the ob- 
served one by Bahcall & Cen (1992a,b). Bartlett & Silk 
(1993) come to a similar conclusion using the Press- 
Schechter (1974) formula. They found that the CDM 
model with the COBE normalization produces a tem- 
perature function of clusters higher than that given by 
the observations by Edge et al. (1990) and by Henry 
& Arnaud (1991); 

— the conflict between the normalization of the spectrum 
of the perturbation which is required by different types 
of observations; in fact, the normalization obtained 
from COBE data (Smoot et al. 1992) on scales of the 
order of 10'^ Mpc requires erg = 0.95 ± 0.2, where trg 
is the rms value of ^ in a sphere of 8/i~^Mpc. Nor- 
malization on scales 10-h50Mpc obtained from QDOT 
(Kaiser et al. 1991) and POTENT (Dekel et al. 1992) 
requires that erg is in the 0.7-^ 1.1, range which is com- 
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patiblc with COBE normalization while the observa- 
tions of the pairwise velocity dispersion of galaxies on 
scales r < 3 Mpc seem to require ug < 0.5. 
— the incorrect scale dependence of the galaxy correla- 
tion function, ^(r), on scales 10-^ 100 /i~^Mpc, having 
^(r) too little power on the large scales compared to the 
power on smaller scales (Maddox et al. 1990; Saunders 
et al. 1991; Lahav et al. 1989; Peacock 1991; Peacock 
& Nicholson 1991); the APM survey (Maddox et al. 
1990), giving the galaxy angular correlation function, 
the 1.2 Jy IRAS power spectrum, the QDOT survey 
(Saunders et al. 1991), X-ray observations (Lahav et 
al. 1989) and radio observations (Peacock 1991; Pea- 
cock & Nicholson 1991) agree with the quoted conclu- 
sion. As shown in studies of galaxy clustering on large 
scales (Maddox et al. 1990; Efstathiou et al. 1990b; 
Saunders et al. 1991) the measured rms fluctuations 
within spheres of radius 20/i^^Mpc have value 2 3 
times larger than that predicted by the CDM model. 

Several alternative models have been proposed in order to 
solve the quoted problems (Peebles 1984; Shafi & Stecker 
1984; Valdarnini & Bonometto 1985; Bond et al. 1988; 
Schaefer et al. 1989; Holtzman 1989; Efstathiou et al. 
1990a; Turner 1991; Schaefer 1991; White et al. 1993a; 
Shaefer & Shafi 1993; Holtzman & Primack 1993; Bower 
et al. 1993). Most of them propose in some way a modifi- 
cation of the primeval spectrum of perturbations. In two 
previous papers (Del Popolo & Gambera 1998a; 1999) we 
showed how, starting from a CDM spectrum and taking 
into account non-radial motions, at least the problem of 
the clustering of clusters of galaxies (first point above) 
and the problem of the X-ray temperature (second prob- 
lem quoted above) can be considerably reduced. 
In this paper we extend the model to two other tests of the 
abundance of clusters: the mass function and the VDF, 
permitting to estimate the expected number density of 
clusters within a given range in mass and velocity, respec- 
tively. We also study the effect of non-radial motions on 
the shapes of galaxy clusters. In recent papers (de Theije 
et al. 1995; de Theije et al. 1997) it has been shown that 
clusters of galaxies are more nearly spherical and more 
centrally condensed than the predictions of CDM models 
with = 1. 

As we shall see, non-radial motions have the effect to pro- 
duce more spherical clusters and are able to reconcile the 
CDM with f2 = 1 predictions on clusters elongations with 
observations. 

In Sect. 2 we shall use the same model introduced by Del 
Popolo & Gambera (1998a,b; 1999) to take into account 
non-radial motions, arising from the tidal interaction of 
the protoclusters with the neighbouring protostructiues, 
and we shall compare the mass function calculated using 
the CDM model, taking into account non-radial motions, 
with the observed mass function obtained by Bahcall & 
Gen (1992) and Girardi et al. (1998). In Sect. 3 we repeat 



the calculation for the VDF and compare the theoretical 
VDF with the GfA data by Zabludoff et al. (1993) and 
with the data by Mazure et al. (1996) and Fadda et al. 
(1996). In Sect. 4 we study the effect of non-radial mo- 
tions on the ellipticity of clusters and finally in Sect. 5 we 
give our conclusions. 



2. Non-radial motions and the mass function 

One of the most important constraints that a model for 
large-scale structure must overcome is that of predict- 
ing the correct number density of clusters of galaxies. 
This constraint is crucial for several reasons. According 
to the gravitational instability scenario, galaxies and clus- 
ters form where the density contrast, d, is large enough so 
that the surrounding matter can separate from the gen- 
eral expansion and collapse. Consequently the abundance 
of collapsed objects depends on the amplitude of the den- 
sity perturbations. In the CDM model these latter follow a 
Gaussian probability distribution and their amplitude on 
a scale R is defined by cr{R), the r.m.s. value of 6, which is 
related to the power spectrum, P{k). In hierarchical mod- 
els of structure formation, like CDM, (j{R) decreases with 
increasing scale, R, and consequently the density contrast 
required to form large objects, like clusters of galaxies, 
rarely occurs. The present abundance of clusters is then 
extremely sensitive to a small change in the spectrum, 
P{k). Moreover, the rate of clusters evolution is strictly 
connected to the density parameter, fio- Then, clusters 
abundance and its evolution are a probe of Oq and P{k) 
and can be used to put some constraints on them. 
The abundance of clusters of galaxies, together with the 
mass distributions in galaxy halos and in rich clusters 
of galaxies, the peculiar motions of galaxies, the spatial 
structure of the microwave background radiation is one of 
the most readily accessible observables which probes the 
mass distribution directly. 

The most accurate way of assessing the cluster abun- 
dance is via numerical simulations. However, there is an 
excellent analytic alternative, Press & Schechter's the- 
ory (Press & Schechter 1974; Bond et al. 1991). Brain- 
erd & Villumsen (1992) studied the CDM halo mass func- 
tion using a hierarchical particle mesh code. From this 
last work it results that the Press & Schechter formula 
fits the results of the simulation up to a mass of 10 
times the characteristic Icr fluctuation mass, M*, being 
M, ~ lOi^6"'^/("'+^)Af0, where b is the bias parame- 
ter and ni is the local slope of the power spectrum. For 
the case of critical-density universes, N-body simulations 
(Lacey & Cole 1994) have been shown to be in extremely 
well agreement with Press & Schecter's theory. 
Press-Schechter's theory states that the fraction of mass 
in gravitationally bound systems larger than a mass, M, is 
given by the fraction of space in which the linearly evolved 
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density field, smoothed on the mass scale M, exceeds a 
threshold Sr'- 



F(> M) = l-erfc ( ^ 

' 2 ' [V2aiRf, 



(1) 



The fraction of the mass density in non-linear objects of 
mass M to M+d M is given by differentiating Eq. with 
respect to mass: 
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n(M, z)dM = - ( - ) ^ ~ cxp(-^)— dA/ (2) 



where pb is the comoving background density, i?f is the co- 

/ \ 1/3 

moving linear scale associated with M, Rf — 
Press-Schechter's result predicts that only half of the mass 
of the Universe ends up in virialized objects but in partic- 
ular cases this problem can be solved (Peacock & Heavens 
1990; Cole 1991; Blanchard et al. 1992). 
The mass variance present in Eq. (|l|) can be obtained once 
a spectrum, P{k), is fixed: 
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where W{kR) is a top-hat smoothing function: 
3 



W{kR) 



{kRY 



(sin kR — kR cos kR) 



(3) 
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and the power spectrum P{k) — Ak"T'^{k) is fixed giving 
the transfer function T{k) : 



T(fc) = Mi±HVM.[l + 1.2fci/^_27fc + 
18 Y a 

+ 347(1 - Va/5)fc3/^ - 18(1 - 0.32a^)k^Y 



(5) 



(Klypin et al. 1993), where A is the normalization con- 
stant, a=(l + z)^^is the expansion parameter and k is 
the wave-number measured in units of Mpc~^. This spec- 
trum is valid for k < 30Mpc^^ and z < 25. The accu- 
racy (maximum deviation) of the spectrum is 5%. It is 
more accurate than Holtzman's (1989) spectrum, used by 
Jing et al. (1994) and Bartlett & Silk (1993) to calculate 
the mass function and the X-ray temperature function 
of clusters, respectively. The spectrum is lower by 20% 
on cluster mass scales than Holtzman's (1989). The spec- 
trum was normalized to the COBE quadrupole Q2 — 17/i 
K, corresponding to trg = 0.66. As shown by Bartlett & 
Silk (1993) the X-ray distribution function, obtained us- 
ing a standard CDM spectrum, over-produces the clusters 
abundances data obtained by Henry & Arnaud (1991) and 
by Edge et al. (1990). This has lead some authors (White 
et al. 1993b) to cite the cluster abundance as one of the 
strongest pieces of evidence against the standard CDM 
model when the model is normalized so as to reproduce the 
microwave background anisotropics as seen by the COBE 
satellite (Bennett et al. 1996; Banday et al. 1996; Gorsky 



et al. 1996; Hinshaw et al. 1996). 

The discrepancy can be reduced taking into account the 
non-radial motions that originate when a cluster reaches 
the non- linear regime as follows. A fundamental role in 
Press-Schechter's theory is played by the value of i5c. This 
value is quite dependent on the choice of smoothing win- 
dow used to obtain the dispersion (Lacey & Cole 1994). 
Using a top- hat window function i5c = 1.7 ± 0.1 while for 
a Gaussian window the threshold is significantly lower. 
In a non-spherical context the situation is more compli- 
cated. Considering the collapse along all the three axes the 
threshold is higher, whereas the collapse along the first 
axis (pancake formation) or the first two axes (filament 
formation) corresponds to a lower threshold (Monaco 
1995). The threshold, 6c, does not depend on the back- 
ground cosmology. 

As shown by Del Popolo & Gambera (1998a; 1999), if non- 
radial motions are taken into account, the threshold Sc is 
not constant but is function of mass, M (Del Popolo & 
Gambera 1998a; 1999): 



1 



rtaL^ • dr 



(6) 



where 6co = 1-68 is the critical threshold for a spherical 
model, Tj is the initial radius, rta is the turn-around ra- 
dius and L the angular momentum. In terms of the Hubble 
constant, Hq, the density parameter at current epoch, f^o, 
the expansion parameter a and the mean fractional den- 
sity excess inside a shell of a given radius, S Eq. (||) can 
be written as (Del Popolo & Gambera 1998a; 1999): 
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where a is the expansion parameter, and a^in its value cor- 
responding to r;. The mass dependence of the threshold 
parameter, 6c{v), was obtained in the same way sketched 
in Del Popolo & Gambera (1999): we calculated the bind- 
ing radius, rb, of the shell using Hoffmann & Shaham's 
criterion (1985): 



Tcir, v) < to 



(8) 



where Tc{r, v) is the calculated time of collapse of a shell 
and ta is the Hubble time. We found a relation between 
V and M through the equation M — A-Kp\,r'^/?>. We so 
obtained 5c[v{M)]. We obtained the total specific angular 
momentum, /i(r, v) — L{r, v)/Msh, acquired during expan- 
sion, in the same way sketched in Del Popolo & Gambera 
(1998a; 1999). Taking account that 5c depends on M Eq. 
(0) becomes: 



n(M, z)dM 



1/2 



p exp(' 



2a^ 



M 



d5_ 



do- 

= dM 



dAf(9) 



The result of the calculation is shown in Fig. 1. Here, the 
mass function of clusters, derived using a CDM model 
with ^Iq — 1, h — \/2 normalized to Qcobe = K 
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Fig. 1. Cumulative mass function calculated using a CDM 
model without taking into account non-radial motions 
(dashed line) and taking account non-radial motions (solid 
line) compared with Bahcall & Cen (1992) data (full dots) 
and with that of Girardi et al. (1998) (open squares). 

and taking into account non-radial motions (solid line), 
is compared with the statistical data by Bahcall & Cen 
(1992) (fuU dots) and with that of Girardi et al. (1998) 
(open squares) and with a pure CDM model with l^o = 1, 
h = 1/2 (dashed line). Bahcall & Cen (1992) estimated 
the cluster mass function using optical data (richness, ve- 
locities, luminosity function of galaxies in clusters) as well 
as X-ray data (temperature distribution function of clus- 
ters). Groups poorer than Abell clusters have also been 
included thus extending the mass function to lower masses 
than the richer Abell clusters. Girardi et al. (1998) data 
are obtained from a sample of 152 nearby (z < 0.15) Abell- 
ACO clusters. As shown, the CDM model that does not 
take account of the non-radial motions over-produces the 
clusters abundance. The result is in agreement with the 
study of the mass function in the SCDM model by Jing & 
Fang (1994) and by Bahcall & Cen (1992a,b). Even with 
a lower normalization, CDM cannot reproduce the cluster 
abundance as stressed by Bartlett & Silk (1993), on the 
contrary a reduction of normalization produces a too steep 
mass function (Bahcall & Cen 1992b; Bartlett 1997). The 
introduction of non-radial motions (solid line) reduces re- 
markably the abundance of clusters with the result that 
the model predictions are in good agreement with the ob- 
servational data. This result confirms what found in Del 
Popolo & Gambera (1999) showing how a mass dependent 
threshold, Sc{M), (dependence caused by the developing 
of non-radial motions) can solve several of SCDM discrep- 
ancies with observations. 



3. Non-radial motions and the velocity dispersion 
function 

The VDF is defined in a similar way to the mass function, 
namely it is the number of objects per unit volume with 
velocity dispersion larger than a^. Since the velocity dis- 
persion CTv can be observed directly (on the contrary, mass 
measurement is usually model dependent), VDF provides 
a good test of theoretical models. Observed CTv comes from 
the measurement of galaxy redshift. The VDF can be cal- 
culated starting from the mass function: 



n(crv) = n{M) 



dM 

dCTv 



(10) 



The cumulative VDF can be obtained integrating Eq. (|10|) : 
n{> (Tv) — / n{av')dav' (11) 



In order to use Eq. (|T^) to calculate the VDF we need a 
relation between the velocity dispersion, CTv, and mass, M. 
This can be obtained in several ways. The typical virial 
temperature may be written as: 



aG^niYi Mvir 

3 -Rvir 



(12) 



where a is a factor of order unity. Evrard (1989) found in 
N-body simulations of a CDM model that typical clusters 
had a ~ 0.8. The molecular weight, fi, corresponding to a 
fully ionized gas with primordial abundances is ~ 0.6 and 
toh is the proton mass. Afvir and i?vir are respectively the 
virial mass and virial radius, which are connected by: 



178— pb^vii 



(13) 



Eq. (|T2|) has been tested in several N-body simulations. 
These numerical simulations give (Evrard 1996): 



Tvh. = (6.8/i2/3keV)A/2/3 



(14) 



The relation is so good that Evrard (1997) uses it as a 
primary mass indicator for observed clusters when con- 
sidering the baryon fraction over an ensemble of clusters. 
Using: 



(15) 



(see Thomas & Couchman 1992; Evrard 1990) togheter 
with Eq. (|l|) we find the necessary relation between cr„ 
and M: 



= 824km/s 



hM 



lOiSMfT 



1/3 



(16) 



(Evrard 1989; Lilje 1990). In Fig. 2 we compare the VDF 
obtained from a CDM model taking account of non-radial 
motions (solid line) with the CfA data by Zabludoff et al. 
(1993) (fuU dots) based on their survey of i? > 1 Abell 
clusters within z < 0.05 and with the data by Mazure et 
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Fig. 2. Cumulative VDF calculated using a CDM model 
without taking into account non-radial motions (dashed 
line) and taking into account non-radial motions (solid 
line) compared with ZabludofFet al. (1993) data (full dots) 
and with those by Mazure et al.(1996) (full triangles) for 
R > 1 clusters and Fadda et al. (1996) (open squares) 
for R > —1 clusters. The theoretical curves are obtained 
using a a^-M relation with zero scatter. 



al. (1996) (full triangles) and Fadda et al. (1996) (open 
squares). Mazure et al. (1996) data are obtained from a 
volume-limited sample of 128 Raco ^ 1 clusters while 
that of Fadda et al. (1996) are obtained from a sample 
of 172 nearby galaxy clusters (z < 0.15). We also plot the 
VDF obtained from a CDM model without non-radial mo- 
tions (dashed line) . The SCDM model predicts more clus- 
ters than the CfA observation except at Cv ~ llOOkm/s. 
As reported by Jing & Fang (1994) the SCDM model 
can be rejected at a very high confidence level (> 6(t). 
The reduction of the normalization reduces the forma- 
tion of clusters, thus resolving the problem of too many 
clusters, but leads to a deficit at ~ llOOkm/s. When 
non-radial motions are taken into account (solid line) we 
obtain a good agreement between theoretical predictions 
and observations. Both CDM and CDM with non-radial 
motions predict more clusters of low velocity dispersion 
{ay < 300km/s) than the observation. This discrepancy is 
not significant because the data at Cv < 300km/s could 
be seriuosly underestimated because: 

— as Zabludoff et al. (1993) stressed, their calculations 
of group velocity dispersions < 300km/s are often un- 
derestimates; 

— Zabludoff et al. (1993) measure only for CfA groups 
with > 5 group members while for CTv < 300km/ s a 
fraction of groups could have less than 5 members. 



Also Fadda et al. (1996) cannot draw firm conclusions 
about their incompleteness level, and hence about the 
behaviour of the a distribution for a < 650km/s, while 
Mazure et al. (1996) has a supposed completeness limit of 
a ~ 800km/s. A comparison between Eq. ( p!^ ) and N-body 
simulations of clusters of galaxies in a CDM model shows 
that Eq. ( p^ holds with a rms scatter of ~ 10% (Evrard 
1989; Lilje 1990). To take into account the scatter in Eq. 
( p^ and Eq. (|l6|) it should be necessary to convolve n(f7v) 
with a Gaussian having dispersion of 10%-;- 20% (see Lilje 
1990). The result of this operation is that the abundance 
of clusters with high velocity dispersion depends on the 
assumed value of the rms scatter: a larger scatter implies 
a larger n(crv). The result of this effect is the worsening of 
the problem of clusters abundance over-production in the 
CDM model without non-radial motions. 



4. Non-radial motions and the shape of clusters 

Most clusters, like elliptical galaxies, are not spherical and 
their shape is not due to rotation (Rood et al. 1972; Gre- 
gory & Tiflt 1976; Dressier 1981). The perturbations that 
gave rise to the formation of clusters of galaxies are alike 
to have been initially aspherical (Barrow & Silk 1981; Pea- 
cock & Heavens 1985; Bardeen et al. 1986) and aspheric- 
ities are then amplified during gravitational collapse (Lin 
et al. 1965; Icke 1973; Barrow & Silk 1981). The elon- 
gations are probably due to a velocity anisotropy of the 
galaxies (Aarseth & Binney 1978) and according to Bin- 
ney & Silk (1979) and to Salvador-Sole & Solanes (1993) 
the elongation of clusters originates in the tidal distortion 
by neighboring protoclusters. In particular Salvador-Sole 
& Solanes (1993) found that the main distorsion on a clus- 
ter is produced by the nearest neighboring cluster having 
more than 45 galaxies and the same model can explain the 
alignement between neighboring clusters (Binggeli 1982; 
Oort 1983; Rhee & Katgert 1987; PHonis 1993) and that 
between clusters and their first ranked galaxy (Carter & 
Metcalfe 1980; Dressier 1981; Binggeli 1982; Rhee & Kat- 
gert 1987; Tucker & Peterson 1988; van Kampen & Rhee 
1990; Lambas et al. 1990; West 1994). Clusters elonga- 
tions and alignement could be also explained by means of 
Zeldovich's (1978) "pancakes" theory of cluster formation 
but this top-down formation model is probably ruled out 
for several well known reasons (Peebles 1993). 
The observational information on the distribution of clus- 
ters shapes is sometimes conflicting. Rhee et al. (1989) 
found that most clusters are nearly spherical with elliptic- 
ities distribution having a peak at e ~ 0.15 while Carter 
& Metcalfe (1980), Binggeh (1982), Plionis et al. (1991) 
found that clusters are more elongated with the peak of 
the ellipicities distribution at e ~ 0.5. More recently de 
Theije et al. (1995), de Theije et al. (1997) re-analyzed 
the data studied by Rhee et al. (1989) and that by Plionis 
et al. (1991) concluding that: 
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Fig. 3. EUipticity, e, of clusters versus mass, M . The 
dashed and solid lines represent e for a CDM without and 
with non-radial motions, respectively. 



a) richer clusters are intrinsically more spherical than 
poorer ones; 

b) the projected elongations of clusters are consistent 
with a prolate distribution with clusters ellipticity dis- 
tribution having a peak at e ~ 0.4 and extending to 
e = 0.8; 

c) in a JIq = 1 CDM scenario clusters tend to be less 
spherical than those in a fio = 0.2 universe and are 
too elongated with respect to real observed clusters. 

To study the effect of non-radial motions on the shape 
of clusters we shall use a model introduced by Binney & 
Silk (1979). In that paper they showed that tidal inter- 
actions between protoclusters and the neighbouring pro- 
tostructures should yield prolate shapes (before virializa- 
tion) with an axial ratio of protostructures of ~ 0.5, the 
typical value found in clusters. After virialization the pre- 
existing elongation is damped and the axial ratio leads 
to values of about 0.7 0.8, that are higher with respect 
to observations. As observed by Salvador-Sole & Solanes 
(1993) this last discrepancy can be removed taking into 
account that tidal interaction keeps going on after virial- 
ization and that on average the damping of elongations 
due to violent relaxation is eliminated by its growth after 
virialization. Then this growth restores a value of e near 
the one that clusters had before virialization. 
According to the quoted Binney & Silk (1979) model, an 
initially spherical protostructure (e.g. a protocluster) of 
mass M having at distance r{t) from its centre a series of 
similar protostructure of mass M' shall be distorted. In 
order to calculate the distortion we must write the equa- 
tion of motion of a particle at position R relative to the 



centre of M (assumed as origin of coordinates). Suppos- 
ing that the effective perturbing mass is less by a fac- 
tor A = (p — Ph)/ Ph than its true mass M', in the limit 
Irl >> \R\ we have: 



R = 



GM 
~R3~ 



R 



GAM' ( Rr 



R 



(17) 



If the tidal interaction is treated as first order, writing 
R{t) = Ro + Ri with Ro = -^^Ro and Rg = Rm.a{t), 
it is possible to show that the component of Ri parallel 
to Rm is: 



3 ( Tfn Rm) / fn 



Rn 



G{x,xi) 



(18) 



where x — 2a and xi is the value of x at which the pertur- 
bation was switched on. The term in square parentheses 
reduces to 2|i?i„|, if Rm is parallel or antiparallel to Tm, 
and to -|i?m| when they are orthogonal. Then the initially 
(at time given by x = xi) spherical density enhancement 
M becomes a prolate spheroid having ellipticity: 

e - (l-b/a) = ^pG{x,x,)/{l+pG) c ^{py/^G{2,x,)il9) 

where G(2, xi) is defined in the quoted paper (see Eq. 13c) 
and {p^y/^ is given by (Binney & Silk 1979): 



2\l/2 _ 



where 



p^{M, M', r)N{M', r)drdA//' 



1/2 



(20) 



Rn 



i{0 

_ ( 3M \ 



M' 

M 
1/3 



being the time of maximum expansion, = '"(im)i 
and N {M\r)dr is the number of condensations of mass 
M' lying between r and r-\-dr from M . We calculated this 
quantity using the Press-Schechter's theory (see Sect. 2): 



N{M',r) ^ 4:nr^n{M') 



(21) 



To calculate n{M') in the case of CDM without non-radial 
motions we used Eq. (B while for a CDM with non-radial 
motions we used Eq. (O) . With the previous definitions 
{p'^y^'^ is given by: 



(p 



2x1/2 



1 



N{M',r)dM' 



M + M' 



1/2 



(22) 



In Fig. 3 we show the shape of e as a function of mass, M. 
As shown e declines with mass in agreement with the 
above quoted point a): richer clusters are more spheri- 
cal than poorer clusters. The physical reason for this may 
be that regions of higher density turn around earlier from 
Hubble flow than lower density regions (de Theije et al. 
1995; Ryden 1995). The fundamental point in which we 
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are interested is the effect of non-radial motions on e. In 
a CDM model that takes into account non-radial motions 
(solid line) e is smaller than in the simple CDM (dashed 
line). This is in agreement with the point c): non-radial 
motions reduce the elongation of clusters. For a cluster of 
IO^^Mq we get a value of e ~ 0.5, if non-radial motions 
are excluded, while e ~ 0.43 when non-radial motions are 
taken into account. Increasing the mass, as expected, clus- 
ters tend to become more and more spherical. Finally Bin- 
ney & Silk (1979) model predicts that even a spherical 
density enhancement M soon becomes a prolate spheroid 
in agreement with point b) and also with Salvador-Sole & 
Solanes (1993) result. 

5. Conclusions 

In this paper, using the model introduced by Del Popolo & 
Gambera (1998a; 1999), we have studied how non-radial 
motions change the mass function, the VDF and the shape 
of clusters of galaxies. We compared the theoretical mass 
function obtained from the CDM model taking into ac- 
count non-radial motions with the experimental data by 
Bahcall & Cen (1992a,b) and Girardi et al. (1998). The 
VDF, calculated similarly to the mass function, was com- 
pared with the CfA data by Zabloudoff et al. (1993) and 
those of Mazure et al. (1996) and Fadda et al. (1996). Tak- 
ing account of non-radial motions we obtained a nothe- 
worthy reduction of the discrepancies between the CDM 
predicted mass function, the VDF and the observations. 
Non-radial motions are also able to change the shape of 
clusters of galaxies reducing their elongations with respect 
to the prediction of the SCDM model. This last result is 
in agreement with recent studies of the shapes of clusters 
by de Theije et al. (1995, 1997). 
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